CASCADE OF PHASE SHIFTS FOR NONLINEAR 
SCHRODINGER EQUATIONS 



REMI CARLES 



Abstract. We consider a semi-classical nonlinear Schrodinger equation. For 
initial data causing focusing at one point in the linear case, we study a nonlin- 
earity which is super-critical in terms of asymptotic effects near the caustic. We 
prove the existence of infinitely many phase shifts appearing at the approach 
of the critical time. This phenomenon is suggested by a formal computation. 
The rigorous proof shows a quantitatively different asymptotic behavior. We 
explain these aspects, and discuss some problems left open. 



1. Introduction 

We consider the semi-classical limit of the Cauchy problem, for (t, x) E K+ x R": 

(1.1) z£d t u £ + yAu e = /(|u e | 2 )u £ ; uf t=0 = e fc / 2 ao(x)e-^. 

In the linear case / = 0, the quadratic oscillations of the initial data cause focusing 
at the origin at time t = 1 in the limit e — > (see Section l2~TJl . In the nonlinear 
case, the effective nonlinear effects strongly depend on the size of the initial data, 
that is on k. Changing notations, we consider: 

(1.2) ied t u £ + ^-Au e = f (e k \u s \ 2 )u e ; uf t=0 = a {x)e- 1 ^ . 

In n, we justified the general heuristics presented in [H], in the case of (|1.2Jl . for / 
homogeneous of degree a, f(y) = y a . Two notions of criticality exist for k: outside 
the focal point, and near the focal point, where the amplitude of u e is strongly 
modified. We described the sub-critical and critical cases. The aim of the present 
paper is to study a supercritical case. 

Consider the case f(y) — y a , and denote a = ka. If ciq E H 1 (M n ) with |x|ao E 
L 2 (R n ) and a < 2/(n - 2) when n > 3, u e is defined globally in time in H^W 1 ). 
The following distinctions were established in £Q: 





a > na 


a = na 


a > 1 


linear caustic 
linear WKB 


nonlinear caustic 
linear WKB 


a = 1 


linear caustic 
nonlinear WKB 


nonlinear caustic 
nonlinear WKB 



The term "linear WKB" means that outside the caustic, the propagation of u £ 
can be described by a geometrical optics approximation, with only linear effects 
involved at leading order. The term "linear caustic" means that nonlinear effects 
are negligible at leading order when the solution crosses the focal point. In either 
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of the two critical cases, nonlinear phenomena are described (the doubly critical 
case was studied more precisely in [2]); we recall the case "nonlinear caustic, linear 
WKB" in Section l2~^l The case we study in this paper corresponds to na > a > 1: 
super-critical caustic with linear WKB regime. 

A look at conservation laws suggests the existence of new relevant scales. In the 
case discussed so far, the conservations of charge and energy read (see e.g. pj]): 

\\u £ (t)\\ L 2 = \\a \\ L 2, 
h\eV x u^t)\\l 2 + J^\\ u %t)\\%tl = const. = 0(1) ~ \\\xa f L2 . 

I (7+1 £^0 A 

When a > na, the boundedness of u £ and eV x u £ in L 2 implies, along with 
Gagliardo-Nirenberg inequalities, for < e < 1: 

e a \\u £ (t)\\l°tl < e™\\u £ {t)\\%tl. < ||u e (t)||^ (n - 3)ff || e V a « e (t)||25 = 0(1) . 

Thus, linear arguments allow us to recover a control of the nonlinear term in the 
energy. Such a line of reasoning fails when a < na: the control provided by the 
conservation of energy hides stronger nonlinear effects. In the linear case / = 0, 
and ao G <S(R n ), one can check that the following point-wise estimate holds: 

(1-4) \u e (t,x)\< -• 

(e+\t-l\) n/2 

In the four cases of the table, the same estimate holds for the nonlinear solution 
in space dimension one (P3E|); like in the linear case, it is sharp. Thus, the above 
estimate e ncr ||u e (t)||^^ 2 = 0(1) is sharp only near the focal point. We now use 
the a priori estimate £ Q ||w e (t)||^2^+2 = 0(1) given by the conservation of energy 
only for t « 1. Assuming that like in all the cases of the above table, at time t = 1, 
u 6 is described by a concentrating profile, 



U £ (t,x) ~ jttS ( — \ , 



we check that the "linear" value 7 = 1 is forbidden (the power of e in front of 4> is to 
ensure the L 2 -norm conservation) . Guessing that the nonlinear term s a 1 1 u e (t) \ \ 2 ^t+2 
in the energy is exactly of order 0(1) at the caustic, we find 7 = a/[na), that is: 

(1.5) 7 =-<l. 

n 

We will not prove that the above argument is correct (see Sectional, but we will 
show that the scale e 7 is an important feature of this problem. Notice also that the 
above argument suggests that the amplification of the solution u 6 as time goes to 
1 is less important than in the linear case; super-critical phenomena may occur in 
the phase, and also affect the amplitude. 

We now go back to the notation 1)1.2(1 . and do not assume in general that the 
nonlinearity is homogeneous (unless it is cubic): 

Assumptions 1.1. The space dimension is n > 2. 

The initial amplitude belongs to the Schwartz space: oq G <S(R n ). 

The nonlinearity is smooth: f G C°°(K+;M). 

/(0) = and f > 0. In particular, the nonlinearity is cubic at the origin. 

Remark 1.2. We suppose ao G 5(M™) to avoid to count derivatives when not neces- 
sary. We could as well assume that ao belongs to Sobolev type spaces. If we require 
a control on the growth of / at infinity, (0 <)/(y) < (y) q for q < when n > 3, 
then for every fixed e > 0, u £ is global in time, continuous with values in H 1 (M n ) 
(see e.g. [H])- This includes a cubic nonlinearity in space dimension two or three. 
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Remark 1.3. The assumption /(0) = is only to simplify notations, since replacing 
/ with / - /(0) turns u £ (t,x) into u £ (t, x)e if( -°W e . 

Remark 1.4. The assumption of the nonlinearity being cubic at the origin is rem- 
iniscent of the paper by E. Grenier |S] (see also P. Gerard 0)- The proof of our 
main result relies on ideas introduced in Sj (see Section [3J. 

Remark 1.5. The one-dimensional cubic nonlinear Schrodinger equation is inte- 
grable. The case k = with more general WKB data was treated in |1()| . 

Before stating our main result, we give the following definition (see e.g. [V2\): 

Definition 1.6. If T > 0, {kj)j>\ is an increasing sequence of real numbers, {fa)j>i 
is a sequence in H°°{W l ) := n s > H s (M. n ), and <j> G C([0, T]; H S (W 1 )) for every 
s > 0, the asymptotic relation 

4>(t,x)~^2 tkj fa ( x ) as 1 -> ° 

means that for every integer J > 1 and every s > 0, 

J 

<^>") _ ^2 tki fa =°{ tkj ) asi^O. 

Theorem 1.7. Lei Assumptions I7"71 &e satisfied. Assume n > k > 1. TTien 
i/iere exist T > independent of e G]0, 1], a sequence (fa)j>i ^ n H°°(W n ) 7 and 
<j> G C([0,T];JJ s (R")) /or every s > 0, suc/i i/iat: 

1. 0(t,x) - Ej>i t^faix) ast^O. 

2. For 1 — t ^> e 1 (7 = fc/n < 1 J, i/ie asymptotic behavior of u £ is given by: 

limsup sup ||u £ (i) - v e (t)\\ L2 — > 0, 

0<t<l-AeT A^+oo 

: 



where v (t, X) = r~7? a ° i ) ex P i£ 



We now comment this result. In the linear case / = 0, the above result holds 
with 7 = 1 and (f> = (see Section ETT) l. We recall in Section 12^1 that in the critical 
case "nonlinear caustic, linear WKB" , the same asymptotic as in the linear case 
holds for 1 — t ^> e. The case k < n is super-critical as far as nonlinear effects near 
t = 1 are concerned. We emphasize two important features in the above result: the 
analysis stops sooner than 1 — t 3> e, and nonlinear effects cause the presence of 
the (nontrivial) phase cf>. For 1 — t 3> e 7 , we have 



1 - t ' 1 - 1 J ^ (1 - t)J' n ~ 1 J V 1 - i 



fe-i 



The above phase shift starts being relevant for 1 — t ~ e"^ 1 (recall that n > fc > 1); 
this is the first boundary layer where nonlinear effects appear at leading order, 
measured by <j>\. We will check that this phase shift is relevant: fa is not zero 
(unless u E = 0, see (|5.1ll below). We then have a countable number of boundary 
layers in time, of size 

1 -t ~ £J"-1 , 

which reach the layer 1— t ~ e 7 in the limit j — > +oo. At each new boundary layer, a 
new phase fa becomes relevant at leading order. In general, none of the fa's is zero: 
see e.g. (|5.3H for fa. The result of a cascade of phases can be compared to the one 
discovered recently by C. Cheverry |S] in the case of fluid dynamics, although the 
phenomenon seems to be different. Yet, our result shares another property with [5], 
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which does not appear in the above statement. Theorem 1 1 . 71 shows perturbations 
of the phase (the <^j's), but not of the amplitude: the main profile is the same as 
in the linear case, that is, a rescaling of ciq. However, to compute the first N phase 
shifts, (4>j)i<j<N , one has to compute N — 1 corrector terms of the main profile clq. 
This appears in Proposition 14. 31 see also Equations (|5.1|) - (|5.3|) . 

The assumption k > 1 means that we start with a linear WKB regime. Indeed, 
for small positive time, u £ remains of order 0(1), and f(e k \u £ \ 2 ) ~ e k \u £ \ 2 f (0). 
The main term is then the same as in pQ with a = 1 and a = k. As recalled in the 
above table, a > 1 corresponds to a propagation which is linear at leading order. 

Each phase shift oscillates at a rate between 0(1) (when it starts being relevant) 
and 0(e 7 ~ 1 ) (when it reaches the layer of size e 7 ). Since 7 > 0, this means that 

each phase shift is rapidly oscillating at the scale of the amplitude, but oscillating 

11 2 

strictly more slowly than the geometric phase 2 J(t-i) ' ^ or 1 — ^ 3> £ 7 - We will see 
in Section that for 1 — t = 0(e 7 ), all the terms in <p, plus the geometric phase, 
have the same order: all these phases become comparable, see 1)6.3(1 . 

We will prove a more precise asymptotics than the L 2 estimate of Theorem II. 71 
see Proposition 14.11 and ((4.1(1 . We restricted our attention to the L 2 norm for the 
sake of brevity. 

Unfortunately, our analysis stops at the boundary layer of size e 7 : we can only 
go up to 1 — t = e 7 /T where T appears in Theorem ll.7l We will discuss this fact in 
Section and explain why we took care of never speaking of "focal point" in the 
super-critical case, but only of caustic (as a matter of fact, even the existence of a 
caustic is not clear, see Section |SJ) . For instance, the geometry of the propagation 
is not known for 1 — t < e 7 , while the analysis shows that it occurs on the rays 
of linear geometric optics before this layer (see Figure ^) . On the other hand, we 
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Cascade of phase shifts 
t along the same rays 



Free transport along rays 
of geometric optics 
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Figure 1. Geometry of the propagation, with morally A — > +00. 



know that the order of magnitude of the amplitude changes to reach the boundary 
layer of size e 1 . Recall that Theorem 11.71 describes the asymptotic behavior of u e 
for 1 — t > Ae 7 , in the limit A — > +00. In this region, leading order nonlinear 
effects are visible only in the phase. As mentioned above, our analysis is valid for 
1 — t > e 1 /T . Between the initial time and this region, the amplitude of u e varies 
like (1 — <) - "/ 2 , and changes from 0(1) to &{e~ k / 2 ) (see also Remark [4.71 for the 
asymptotic behavior of u e for t < 1 — e 7 /T). 

The rest of the paper is organized as follows. In Sectional we recall the proof of 
the analog of Thcorcm ll.7l in the linear and critical nonlinear cases. In Section |3 
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we present a formal computation that suggests a result like Theorem II. 71 Based on 
the result by E. Grenier [H] and a "semi-classical conformal transform" , we give the 
proof of Theorem 1 1.71 in Section 0] In Sectional we compare the rigorous approach 
with the formal result of Sectional The discussion about some possible phenomena 
for t > 1 — e 7 /T appears in the final Section 



Acknowledgments. The author is grateful to Christophe Cheverry for stimu- 
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2. Free and critical cases 
2.1. The linear equation. Consider the linear equation: 

(2.1) ied t uf in + ^Auf in = , (t,x)6lxl»; < l|(=0 = ao^e^ . 

As e — > 0, the rays of geometric optics (classical trajectories) are lines = const., 

and meet at the origin at time t = 1. Indeed, the bicharacteristic curves are defined 

If I 2 

by the Hamilton flow associated to p(t, x, r, £) = r + 4j-: 

i=\ ■ x = £ ; r = C = ; ar(0) = x ; f(0) = V0(O,a(O)) = -x ■ 
Of course, uf in can be expressed in terms of an oscillatory integral: 

(2.2) u^x)^^^l^^a (y) d y. 

Applying stationary phase formula yields the same result as using WKB methods 
below, up to the same boundary layer. Seek 



uf in (t,x) ~ vf in (t,x) = v°(t,x)e t 

e—>0 

Plugging this into 1|2.1|1 and canceling the 0(e°) and 0(e 1 ) terms, we find: 

9^+~|V^| a =0, <£(0,z) = -^; 

d t v° + V x <p ■ V x v a + iu°A0 = , v°(0, x) = a (x) . 

For t < 1, one has explicitly: 

\x\ 2 

^>(*,a;) = 7777 77 ! = 7" 77T77 a 



2(t-l) ' yi ' (1-t)"/ 2 \l -t 



Moreover, vf in solves: 



2 lm 2(1 -i) 2 (1 -t)»/2 V 1 -^ 

£ 

(T-t)2 



Let r £ (t,x) denote the source term: ||r e (£)||£2 < — — -. Standard energy esti- 



mates for Schrodinger equation yield: 



ei\\» e (t)-if(t)\\v<V(t)\\L* < ' 



dt" w y ;iUj ~ (1 -t) 2 ' 

and 

sup \\u s (s)-v s (s)\\ L2 <- f \\r%,s)\\ L2 ds < 

0<s<t £ Jo 1 — I 
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Thus, WKB approximation is interesting up to a boundary layer in time of size e 
before the focus. Moreover, vf in is exactly the approximate solution of Theorem 1 1.71 
with = 0. Past this boundary layer, l|2.2(l shows that for |1 — t\ = 0(e), 

■ i ® i 2 ■ i x i 2 

where J 7 denotes the Fourier transform. For ( - 1 > e, stationary phase formula 
yields the same asymptotic description as above, up to the Maslov index (see |Hin])- 
In particular, we see that (|1.4|l holds for uf in , and is sharp. 

2.2. The critical nonlinear case. We recall the main result of . Consider (|1.2() 
in the case k = n, which is critical concerning the role of the nonlinearity near the 
focal point. Introduce the scaling 

Then the function ip e solves 

(2.3) id t r + ^r = f(\r\ 2 )r ; 0L=i = £ n/2 a (ex)e- isl ^ . 

Z I e 

A way to understand criticality is that e has disappeared from the equation satisfied 
by ip E . Using global well-posedness results for nonlinear Schrodinger equations 
(under assumptions on the nonlinearity which are different from Assumptions 11.11 
see e.g. PJ), one has 

W - ^llz,-(R;Hi(R»)) -> as £ ^ , 
where -0 is the (global) solution of the Cauchy problem 

(2.4) zd t 0+iA^/(H 2 )0 ; e- i i A i>(t,x)\ t= _ oo = F-\a )(x). 

Scattering theory shows two interesting features: for large i/j(t,x) behaves like 
a solution of the linear Schrodinger equation. This implies that for 1 - t > £, the 
solution u 6 can be approximated by uf hl (or vf-): no nonlinear effect is relevant 
before the same boundary layer as before. The second point is that for |1 — 1\ < s, 
nonlinear effects occur at leading order, and are measured (in average) by the 
nonlinear scattering operator associated to (|2.4() . 

3. A FORMAL COMPUTATION 

From now on, we assume that k < n. To simplify notations, and since the 
Assumptions 11.11 will be needed for rigorous proofs only, consider the case of an 
homogeneous nonlinearity: f(y) — y a , and denote a = ka. Then i|1.2|) becomes 

(3.1) ied t u e + ^-Au £ =e a \u e \ 2a u e ; uf t=0 = a (x)e~ 1 ^ . 

The caustic is supercritical: no > a. We also assume a > 1 (linear WKB). Because 
this section is only formal, we shall be very brief about the computations, and only 
give the main steps. 
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3.1. A first boundary layer. Two approaches (at least) lead to the same result: 
Lagrangian integral with stationary phase formula (like in where the critical 
one-dimensional cubic case is considered), and generalized WKB methods. We 
shall retain the second one, which we use in the next subsection. Seek 

u £ (t,x) ~ v{{t,x) = u°(t,x)e l 

e— >0 

and change the usual hierarchy to force the contribution of the nonlinear term to 
appear in the transport equation: 

^+i|V,0| 2 = O, 0(0,*) = -^; 

d t u° + X7 x <f) ■ X7 x u a + ^u°Acj) = -i£ a - 1 |u | 2tr u , u°{0,x) = a (x) . 

The eikonal equation is the same as in Section 12.11 as well as its solution. The 
transport equation is an ordinary differential equation along the rays of geometric 
optics J--7 = const., of the form 

y = -ie a - 1 \y\^y. 



The modulus of u° is constant along rays, and 
u°(t,x) = (1 _ 1 t) „ /2 «o exp I -ie a ~ l 



«o 



1 - t 



2 ° " 4 ds 



o (l-s)" c 



Note that the notation is no longer relevant, since u° now depends on e. We have 
a new boundary layer in time, of size before the focus, where 

no — 1 

For 1 — t ~ e@, the above phase shift measures relevant nonlinear effects. We have: 
with 



\\rl(s)\\ L 2ds 



< 



e J »■ xv-ziix. ~ 1 (l-£)2na-l- 

Following the energy estimates of Section |2~T1 this quantity might be the one that 
dictates the size of the error u £ — v\ (see Section [E] for a discussion on that issue). 
The second term is "new" (the first term is the same as in Section ing , and suggests 

2a-l 

the existence of a second boundary layer, of size e 2 ""- 1 . 

3.2. Infinitely many boundary layers: cascade of phase shifts. Seek an 
approximate solution of the form: 

We find 

2/2 2 

ied t v e + £ —&v £ = !yA s £ - ed t g £ - £ —\V x g E \ 2 + j^-jx ■ V x g* 

+i (l4fTT v ^ ' Va ° (i^) ^ + \ (l^) 2 (T^ Aao (j4t 

As suggested by the previous paragraph, write 

(3.2) g e (t, a) = ^ * ( (1 _^ )ng 'T^) ds » with ~£^<&(0 
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In the equation solved by v e , the last term is the "same" as in the linear case: it 
becomes relevant only in a boundary layer of size e. Since our approach will lead 
us to the boundary layer of size e 7 (recall that 7 = k/n = a/na < 1), we ignore 
that term. 

The remaining terms with a factor i are of order, in L 2 , 

sWWII- + r^llVx* e (*)IU- * { (I^^ (T$^' 

and their contribution is also left out in this computation. 
Now we require that v £ be an approximate solution to (|3.1|l : 

(3.3) [d t --^—.V x )g s + ~\V x g e \ z = -- ' 



«o 



1 - < 



i-t y a 2 (1 -fr^ 

Using (|3.2|l . we get: 

5l (£) = -M£)| 2CT , 

(3 - 4) fOT j > 2 ' "® - 4 JC (^-dU-D ^ ' Vff * ' 

with the convention go = 0. This algorithm produces smooth solutions provided 
that |ao(£)| 2CT is smooth (a £ N* or do Gaussian for instance). We neglected the 
terms corresponding to s = in the integration 1)3. 2[) : this does not increase the 
error, since na > a > 1. Defining 

N /"* / £ 7 \ ™ CTJ 

<?w(*> a;) = 7 I] / 1 1 ds x 3j 



s / V 1 - t 



the approximate solution ^ solves 



e 



2 



ied t v% + — Av £ N = s a \v e N \ 2a v e N +7% ; v £ N{t=0 = a (x)e 



2 

with, for 1 — t > e 1 : 



I ft £ (N+l)a~l 

'\r e N (s)\\ L 2ds < 



£ Jo ~ (1 _ t)(iV+l)n CT -l (l_t)TK7- 

To compare with Theorem 11.71 remove the terms corresponding to s = in the 
integration (recall that no > 1), and define: 



-8 J * J V 1 - t 



By definition, we have — (*) Hi 2 = 0(e Q_1 ) for 1 — t > e 1 . One can check 



that v e N solves 



ied t v% + £ -^Av% = s a \v s N \ 2 °v% + r% , 
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with 



r E N (t,x) = [- — - v E N {t,x)g N 



N+l 



(3.5) 



2\l-tJ (1 - t) n / 2 \1 - t 



\ , / x 3 



(l_f)n/2 W U \1 - t 

We have the following result: 

Proposition 3.1 (Formal approximation to l|3.1|l ). Let na > a > 1, do S 

anci fix N eN*. Denote 

e 3 " -1 1 / a; \ 

^(t, ar) = ^ (1 _ t) in*-i inff _ia ^1^7 J ' 

and let v E N be the associated approximate solution. The function v 1 ^ solves 

isd t v% + ^-Av £ N = e a \v s N \ 2a v s N + r% ; i^ |t=0 = a Q (x) e -^ + (e"" 1 ) m L 2 . 

For 1 — t > e 1 = e™, the source term satisfies: 

I ft £ (N+l)u-l £ a 

'\r e N (s)\\ L 2ds < 



e J Q 11 Ny n]L ~ (1 - t)(W+l)n«x-l 

For 1 < j < N, the j term of the series defining g e N becomes relevant in 
a boundary layer of size e n i°- 1 : in the limit N — ► +oo, a countable family of 
boundary layers appear, between s@ and e 7 . In the case cr = 1, which is the only 
homogeneous nonlinearity consistent with Assumptions 11.11 we have a = k and we 
find the boundary layers announced in the introduction. 

Letting N — > +oo (using Borel lemma, see e.g. we find: 

^V(*)IM*< (r^p, 

which is small for l-i>e 7 . 



Remark 3.2. In the critical case a = rwr > 1, we have f3 = 7 = 1: the above 
boundary layers "collapse" one on another. There are no such phase shifts as 
above. 

We point out that the sole estimate of the source term proves nothing. In a 
stability argument, the nonlinearity |u e | 2cr u' r — \v e N \ 2a v e N is usually treated by a 
Gronwall type argument. If the nonlinearity is "too strong", then the above esti- 
mate, which is completely relevant in the linear case, does not necessarily account 
for the size of the error. Since we are in a super-critical case, it is not surprising 
that Proposition 13. II is only a formal result. This remark can be compared to the 
approach in [Jj. To justify a WKB expansion for the nonlinear equation 

ied t u s + ^-Au s = f(\u s \ 2 )u £ , 
constructing an approximate solution that solves 



£ 2 

ied t ut pp + yA< pp = / (|< pp | 2 ) < pp + (e°°) , 



10 



Ft. CARLES 



is not sufficient. Indeed, the computations in [7] show that energy estimates and 
Gronwall lemma do not yield better than 

h e (t)-< pp (t)\\ L2 <e ct ^0(e°°). 

This is the reason why in 7 , WKB expansions are justified for analytic data. This 
assumption yields a source term for u| pp which is 0(e~ 5 / £ ), counterbalancing the 
exponential growth of Gronwall lemma to lead to a good approximation on [0, T] 
for T > independent of e. 



4. Rigorous results 



We now prove Theorem II .71 We will see that the approximate solution we find 
diverges from the one constructed above, a fact which we discuss in Section [5j 

4.1. Semi-classical conformal transform. Introduce the new unknown function 
■0 e given by: 



(4.1) 



u e (t, x) 



1 



i-t \-t 



(1 - i)«/a 
Recalling that 7 = £ < 1, denote 
(4.2) h = 

Changing the notation V' £ ( T i£) hito ip h (t,x), we check that H1.2(l becomes: 
H 2 



.1-7 



0. 



(4.3) + -AV> n = i- 2 / (nv/f ) ^ ; V R | 



a (x) . 



The singular term t 2 in factor of the nonlinearity is actually harmless: as t goes 
to zero, t- 2 f (t n \ip h \ 2 ) « t"- 2 |V> ft | 2 /'(0), and is bounded since n > 2. 

The proof of Theorem 1 1.71 is now reduced to the asymptotic expansion for ip h as 
h — > for i G [fi- 1 ^, j]- Denote ig = ftJ^i . We shall prove the following: 

Proposition 4.1. Let ^4ssMmptaon.s r7~n be satisfied. Assume n > k > 1, and let 

s G N. Then there exists T > independent of H G]0, 1] swc/i that for t G [£o,T], 
the function ip h can be written as ip n {t,x) = a h (t,x)e l< ^ f*^)/* 8 , with 



a \\L°°([t%,T];H s ) + 



lL°°([tg,T];if«) 



as h — > . 



where (a, 0) solves 
(4.4) 



9*0+ -IV^^I 2 +i- 2 /(*"|a| 2 ) -0 
3 f a + V x ■ V^a + ^-aA(j) — 



'|t=o 



0, 
a . 



Moreover, 



lim sup sup 

fi^O t^<t<T 



^ H (t,x)-a (x) e l ^ t > x)/h 



asr^O. 



The second point of Theorem 11.71 follows from the above proposition, since the 
transform (|4.1|l is L 2 unitary (see Proposition 14. 31 below for the asymptotic expan- 
sion of 4>)- We could also include not only derivatives in the above estimates, but 
also momenta. As announced in the introduction, we chose to leave out this refine- 
ment. Note that except for two aspects, Proposition 14. II is nothing but rewriting 
Theorems 1.1 and 1.3 of In our case, time is present in the nonlinearity, and 
data for ip h are prescribed at time ig (with ig — > as H — > 0) instead of time zero. 
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4.2. Construction of solutions to (|4.3() . We recall the ideas introduced by 
E. Grenier [H], and show how to handle the presence of time in the nonlincarity. The 
main idea in [2] is to write the solution of a semi-classical nonlinear Schrodinger 
equation as a WKB solution, where not only the amplitude may depend on the 
small parameter, but also the phase. This changes the usual WKB hierarchy, and 
overcomes the difficulties pointed out in |7| . Seek the solution of (|4.3|l of the form 



ip H (t,x) = a h (t,x)e i * R ( t > x V h , 



with 



(4.5) 



d t J> h 



n\h\2\ 



- W 1 +t- z f(t n \a 







d t a h + V(j} h ■ Va h + ^a h A(t> h = i^Aa h 



0, 
a . 



t— 1 



Introducing the "velocity" v = V<fi h , (|4.5(1 yields 

d t v h + v h ■ Vv h + 2t n - 2 f (t n \a h \ 2 )ReC^Va h ) = ; v h \ t=th = . 

( 4 - 6 ) in 

d t a h + v h ■ X7a h + -a h divv h = i-Aa h ; a 
2 2 

To force the initial time to be zero, introduce 

v n (t,x) = v h (t + t%,x) ; a h (t,x) =a h (t + t%,x) 
Then 1)4. 6|) becomes 

3 t v h + v h ■ Vv n + 2 (t + t%) n ~ 2 f ((t + 4) n |a ft | 2 ) Re (^VS n ) 



0; v* 



\t=o 



0, 



d f a h 



v h ■ Va h + ^a h divv H = i^Aa h ; a h \ t=Q = a . 



Notice that if n = 2 and /' = const. (2D cubic equation, which is conformally 
invariant), the above system is exactly the same as in [H]. 

Separate real and imaginary parts of a , a h = a fi + ib . Then we have 



(4.7) 



d t u h + J2 A M h )d 3 u h = -Lu h . 
/ a* \ 



with u 



L 




Onxn 



Th 

and A(u,0 = J2 A i( u )tj = 



n-2 



where /' stands for /' ((i + tg) n (|S| 2 + |b| 2 )). The matrix A(u,£) can be sym- 
metrized by 



S = 



h 








which is symmetric and positive since /' > 0. We now reproduce the ideas of [5], 
inspired by hyperbolic theory, see e.g. [TT]. For an integer s > 2 + n/2, we bound 
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(Sd"u h , 9"u") where a is a multi index of length < s, and (•, •) is the usual L 2 
scalar product. We have 

j (Sflgu* <9> n ) = (d t Sd%u H , d*u h ) + 2 (Sd t d%u n , d*u h ) 

since S is symmetric. For the first term, we must consider the lower n x n block in 
S. Differentiating (t + tg) 2- ™ yields a non-positive term, and we get 



(d t Sd a x u\ d%u h ) < jjd t (/' ((* + t h ) n (|£ R | 2 + |bf))) ^ (SdZu\ d«u h ) 

So long as ||u ||x,<» < 2||ao||_L°°, we have, for t < 2 (to fix the ideas), 

/' ((i + (|af + |bf)) > inf {f(y) ;0<y< 2" +2 |Kll!oo } = S n > , 
where S n is now fixed, since /' is continuous with /' > 0. We infer, for t < 2, 



U(/'((t + «) w (Ia»| a + |b«| a ))) 



<C(\\u n \\ L ~>)\\d t u h \\ L ~<\\u h \\ H s, 



where we used Sobolev embeddings and i|4.7[) . For the second term we use 
h - 

(Sd t d2u\d2u h ) = - (SX(9> ft ),e?> fi ) - (s82(£M" h )9i" n ),%* f 

3=1 

We notice that SL is a skew-symmetric second order operator, so the first term is 
zero. The second term can be rewritten under the form 



(5a?(53^(u fc )a J -u h ),a?u fc ) = (sE^( uft )^ 9 > R > a > fi ) 

i=i 

n n 



By symmetry of SAj(u ), 



i=i 

n 



Therefore, so long as ||u r ||loc < 2||ao||L°° for t < 2, 



< ll^u^ll 2 IIV II s II < 1 1 -ii^ 1 1 3 



The usual estimates on commutators (see e.g. |TT]) lead to 



i=i i=i 
Notice that 5 _1 can be bounded by C(||u||h«), thus we have proved: 
d 
dt 



<c(h h \\ H s)h h C« 



- £ (Sd^u H ,d^u n )<C(\\u h \\ He ) (SdZu\d%u n ) , 

\a\<s 



\a I <s 



for s > 2 + d/2. Gronwall lemma along with a continuity argument yield the 
counterpart of Theorem 1.1]: 
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Proposition 4.2. Under Assumptions ^! . l\ with n > k > 0, let s > 2 + n/2. Then 
there exist T > independent of h g]0, 1] and ip h (t, x) = a h (t, x^e 1 ^^ 1 '^^ solution 
to O on [t%,T + t%\. Moreover, a h and (j) h are bounded in L°°{[t%,T + tg]; H s ), 
uniformly in h s]0, 1]. 

4.3. Convergence and small time properties. We can now complete the proof 
of Proposition 14.11 For s > 2 + n/2, we know that a h and v h are bounded in 
L°°([0,T];H s ), uniformly in h G]Q,1]. Using g2J, we infer that d t a h and d t v h 
are bounded in L°°([0, T]; H s ~ 2 ). Therefore, a subsequence of (a fi ,v R ) converges 
uniformly in C([0, T]; Hf oc ) to (a,v) solution of (|4.4Jl for any s' < s — 1. Decreasing 
T if necessary, we know that (|4.4|) has a unique solution in C([0, T];iF) ( see e.g. 
By uniqueness for 14.4|l . the whole sequence (a h ,v h ) is convergent. 
For t% < t < T, write 

v h (t, x) = v h (t, x) - [ d t v h (s, x)ds ; a h {t, x) = a h (t, x) - f d t a h {s, x)ds 

Still using the boundedness of dta h and dtv h in L°°([0, T]; -ff s ~ 2 ), we deduce that 
the sequence (a R , v n ) also converges to (14.4(1 . To estimate the convergence and prove 
that it holds not only in H* oc , but also in H s , we shall use the analysis of 1(4. 4(1 for 
small times. 

Proposition 4.3. Let Assumptions^^ be satisfied, with n > k > 1. Then there 
exist sequences {4>j)j>i an d ( a j)j>i * n H cc (W n ), such that the solution of ((4.4(1 
satisfies 

4>{t, x) ~ tjn ~ 1( t ) j( x ) i anrf ~ tjn<1 o( x ) ast^O. 

j>l j>0 

Plugging such asymptotic series into 1(4.4(1 , a formal computation yields a source 
term which is 0(t°°) as t — > 0. The first terms are computed in 1(5.1(1 - 1(5.3(1 . Note 
that the series for a starts with j = 0: the notations are consistent. The result then 
follows with the same approach as in the proof of Proposition ^. 21 and Borel lemma 
(see e.g. We can now be more precise about the convergence of (a h ,4> h ): 

Proposition 4.4. Let s 6 N. There exists C independent of h such that for every 
4<t< T, 

\\a h (t) - a(t)\\ H s + ||^(t) - #t)|| H . <C'(ht + h 2 ^) . 

Proof. We keep the same notations as in the previous subsection, 1(4.7(1 . Define 
(a, v) from (a, v) by the same shift in time as above. Denote by u the analog of u h 
corresponding to (a, v) . We have 

n n . 

dt (u n -u)+J2 AM h )di (u fi - u) + £ (M* h ) A-(u)) 0,u = 2 iu " ■ 

i=i i=i 

Keeping the symmetrizer S corresponding to u h , we can do similar computations 
to those of the previous paragraph. Note that we know that u h and u are bounded 
in L°°([0, T];H S ). Denoting w R = u ri - u, we get, for s > 2 + n/2: 



j t J2 (Sd°w\d«w h )< J2 {Sd^\d^ h )+h\\w\t)\\ Hs . 

\a\<s 

Using Gronwall lemma and Proposition 14 . 31 we infer: 

l|w ft (i)ll*- <^+|(a ) v)| t=tS -(a,v)| t=0 | ffs <M+(4) n - 1 . 
This completes the proof of Proposition 14.41 □ 
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So far, we have not used the assumption k > 1. It appears when one wants to 
approximate e 1 ^ l h by e % ^l h : the factor 1/h requires some care. Proposition 14.41 
shows that for t% < r < T, 

sup \\<f> h (t) - 4>(t)\\ H s (Rn) = (h 2 ^) + h0(r) . 

tg<t<T V ' 

Recalling that 7 = k/n, we then have, for s > n/2, 



+ h ^IL°°([t»,T];ff») 

< 0(h) + 0(r") + + 0(t) . 

This completes the proof of Proposition 14.11 and Theorem 1 1 . 71 follows . 

Remark 4.5 (Well-prepared data). If we had t[} = 0, then the assumption fe > 1 
could be weakened to k > 0. Back to the transform 1)4. 1(1 . if we assume that 

— • i x i 2 _i 
uf t=0 = ao(a;)e 1 2e exp (ie 7 </>(e 7 ,:r)j , 



fc-i 



then the ^ 1_ ' Y J term in the above estimate disappears, and we can conclude 

as before, supposing only n > k > (but still n > 2). Recall that if < k < 1, 
then nonlinear effects are relevant at leading order for any positive time (nonlinear 
propagation); they show up precisely in the phase <f>. 

Remark 4.6. The cascade of phase shifts can be understood as the creation of 
a new phase, appearing discretely in time. With the transform 1(4. 1[) in mind, 
the asymptotic expansion of the phase shift <fi stems from Proposition 14.31 The 
coupling in 14.4(1 shows that even if <j)\t=o — 0; 4>(S, x) is not identically zero, for 
5 > arbitrarily small. The phase of ip h is given asymptotically by 



11-1 



With the same line of reasoning as in the introduction, a phase shift appears for 
t of order h 1 /^- 71 ^ 1 ^, then a second for t of order S 1 ^ 2 ™ -1 ), and so on. The super- 
position of these phase shifts, which are oscillating faster and faster, finally leads 
to a continuous phase, corresponding to an oscillation associated to the wavelength 
h. The idea that the cascade of phase shifts corresponds to the (discrete in time) 
creation of a new (continuous) phase is reinforced by Equation 1(6.3(1 below. 

Remark 4.7. Like in [Hj, we can prove a WKB type asymptotics, not only for small 
times as in Proposition 14.11 but on the whole interval [0, T]. Indeed, with the 
notations of the proof of Proposition 14.41 w fi = 0(h) in L°°([0, T]; H s ) for any 
s £ N. Reasoning as above, w h /h is bounded in L°°([0, T}; H s ), hence dtw h /h is 
bounded in L°°([0, T}; H s ~ 2 ), and a subsequence converges to the linearization of 
(|4.7(l about u. By uniqueness, the whole sequence is convergent, and we have an 
error estimate in the same spirit as in Proposition ^. 41 Denoting d = lim ^ ^ , we 
have: 

ip - ae iA e^/ n — ► . 

L°°([0,T];H s ) fnO 

This is a WKB type asymptotics, with the amplitude ae ld and the phase <fi. 
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5. Stability issues 

The construction of Section |31 and the results of the previous paragraph do not 
agree. To see this, we come back to Proposition 14.31 in (|4.4(l . we have 

(5.1) &{t n - 2 ): Mx ) = -±-f'(0)\a Q (x)\ 2 , 

(5.2) 0(t n - 1 ): ai + V^i • Va + ia o A0! = , 

(5.3) 0(i 2 "- 2 ): (2n-l)0 2 + i|V0 1 | 2 + 2Re(^a 1 )/'(O) + ^|a o | 4 = O. 

The function <f>i is the same as the one obtained by the approach of Section |3J the 
two approximate solutions are close to each other up to the first boundary layer, 
when the first phase shift appears. On the other hand, we see that to get fo, 
the modulation of the amplitude (ai) must be taken into account; in 1)3. 3)) . gi is 
computed without evaluating Aao, unlike (f>2- This means in particular that the 
two approximate solutions diverge when reaching the second boundary layer: the 
approach of Section [3] is only formal, and does not lead to a good approximation. 
And yet, the source term in Proposition 13.11 is small: thus, the linearized semi- 
classical Schrodinger operator is not stable, in the semi-classical limit. We will see 
below that this instability is not due to a spectral instability, but to the fact that 
the approach followed to construct the formal approximation was too crude. 

This phenomenon is due to the super-criticality of the problem. Indeed, for 
fixed e, we deal with a nonlinear Schrodinger equation with repulsive nonlinearity 
(/' > 0), for which global well-posedness results are available (see Remark 11.211 . 
When using the transform (|4. 1|) , notice that the parameter h in 1)4. 3J) goes to zero as 
e — > only when n > fc, that is in the super-critical case (compare with Section ^. 2|l . 

To understand better the instability mechanism, let us go back to the comparison 
between the construction of Section [5] and the results of the previous paragraph. 
Letting iV — > +oo in Proposition ^. II we have an approximate solution of the form 

g l 2e(t-l) / x \ ( 1 — t 

w£ (*> x ) = 75 Iw? a o i 1 ex P l - 



(1-t)"/ 2 Vl-V V e "\l-tl-t 
7 a Q (0 exp - 



T 



(^Mf^t'T^)' 



This formula and the transform (|4.1|l show that the approximation of Section |21 
is too crude, since it ignores the coupling between phase and amplitude for 1)4. 3fl . 
Proposition ^. 31 and (|4.4|l show that to have a good approximation of the phase, the 
coupling between phase and amplitude must be taken into account at every order. 

We can go one step further in the understanding of this apparent instability, 
by applying the transform (|4.1J) to the intermediary approximate solution v e N . We 
show that the formal approximation stops being a good approximation between 
the first and the second boundary layer. Assume a = 1 so that the homogeneous 
nonlinearity satisfies Assumptions 11.11 Like for the exact solution, write 

v 6 N (t,x) = -—i^^N ' e ^(t-D . 

NK ' ' (l-i)"/ 2VJV \l-t l-tj 

Using the expression (|3.5)l . we check that ip% solves 

mi>% + y = t n ~ 2 \iP% | Vat + 4 (*, x) , 
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along with the initial condition tp%\ _ ^_ = a (x) + (ft( Q - 1 )( 1 ~T)) in H s (R n ) 
for any s > 0, where: 



1n\ 



T (t, x) = (t^ N+1 ^ n - 2 K Q {x) + ihKi(t, z)) ip%(t, x) + ihK 2 (t, x) + h 2 K 3 (t, x) , 
for some "nice" functions Kj. Now write Jp%(t,x) = a r ^(t,x)e 1 ^« < - t ' x ^ h . We have: 



d t w h + Mv h )d 3 v h = -Lv h + S h (t, x) , with v h (t, x) 



di<t>% 



3 = 1 



\ d n <f>% J 



(5.4) 



and S n (t, x) = (t + t%) 



( K x + Re ((K 2 - ihK 3 )e l4, »/ h ) \ 
K x + Im UK 2 - iHK 3 )e^/ h 



J 



where the matrices Aj are the same as in Section 14.21 and the functions in the 
definitions of v h and are evaluated at (t + t$,x). We can proceed like in Sec- 
tion l4.2l the new term is the source S r \ Unlike for the exact solution, the oscillatory 
aspect of the problem has not disappeared: the first two components of S h con- 
tain a highly oscillatory factor. Therefore, we cannot expect H independent energy 
estimates here. To measure the effect of this oscillatory term, forget the shift in 
time, and take <g = 0. Then assuming that for small times, dx4>%(t,x) = 0(i n_1 ) 
for any multi- index a (like for the exact solution), the H s norms of the first two 
components of S'' are controlled by 

j.l+s(n— 
0[t+--—- 



A source of order & {t) is not a problem, since we eventually consider the limit t — > 0. 
On the other hand, let us examine the last term. Back to the initial variables, this 
yields a control by 



l+s(n-l) 



e 



-s(l- 7 ) _ 



\-t / 

This is small for 1 - t » e u , with 

7 + sa 

LJ = ; 

1 + s(n- 

We check that for n > a = fc > 1, we have 



(1 - t)i+«(»-i) 



1) 



a a ~ 1 7 + s ( a - 1) 2a - 1 

P = 7 <^ = T— - f 7T < 7, 7 . for any s > 1 . 

7i — l 1 + s(n — 1) 2n — 1 

The first inequality means that we can expect the formal approximation to be a good 
approximation of the exact solution beyond the first boundary layer (and indeed, 
it is close to the approximate solution of Section |2J|. The second one explains why 
the approximation ceases to be relevant before the second boundary layer. 

A possible way to understand the above computation is that the choice of the 
variables is crucial: working with the "usual" unknown v £ (as in Section is not 
very efficient. On the other hand, with the variables introduced by E. Grenier for 
his generalized WKB methods, a precise and rigorous analysis is possible, via the 
transform (|4.1() . Thus, adding new variables helps the analysis: this goes in the 



CASCADE OF PHASE SHIFTS FOR NLS 



17 



same direction as the general theory of geometric optics, and the recent approach 
followed by C. Cheverry for a refinement of this principle 

6. After the cascade of phase shifts 

As announced in the introduction, our analysis stops for times of order t = 
1 — Ae 7 . For A — * +oo, we have Theorem 11.71 For bounded A, A £ [y,+oo[, 
the first part of Proposition 14.11 provides an asymptotic description; see also Re- 
mark l4~7l This shows in particular that the solution tp h is approximated in terms of 
nonlinear geometric optics: the eikonal equation contains the amplitude, therefore 
the geometry of propagation needs not be the same as before, which occurred along 
ra y s tt* = cons t- Note also that the transform 14.1fl changes the space variable 
into a parameterization of the family of rays, when they are straight lines. This 
explains Figure ^ Moreover, between the initial time t = and t = 1 — the 
order of magnitude of u e changes. Indeed, Proposition ^, ll shows that for t £ [£q, T], 
i/j h is of order 0(1) in L°°(R ra ) (take s > n/2 and use Sobolev embeddings). By 
(14.1(1 . we infer that the amplitude of u e varies like (1 — t)~ n / 2 , and changes from 
0(1) initially, to 0{e~ k / 2 ) for 1 — t w e 7 . Such an amplification is similar to what 
happens in the linear case. 

The semi-classical conformal transform (I4.1jl cannot be interesting for values of 
t too close to 1, since it becomes singular. It seems reasonable to introduce the (L 2 
unitary) scaling transform, 

,r, N r, S 1 r ft~ 1 X \ 1 / 1 — 1 X \ 

6.1 u e {t,x) = jztp? ,_ . 

y ' y ' e n ~fl 2 r \ £ 7 e~< J e k / 2 V £l £l J 

With the same change of notation as for ip in Section l4~Tl we have 
(6.2) iW tV h +^A^ = f{\^\ 2 )^. 

We now have exactly the same equation as in [H] . On the other hand, let us examine 
the initial condition. Taking into account the data u e (0,x) = ao(x), we find 

that one may try to decouple to 

I / X \ I a? 1 ^ 

ip h [t, x) ~ , , ,„ a ( - ) e 1 ^ 7 " as t — > -oo and h — * . 
Y v ; | t |n/2 u \ t ) 

Of course, the above limits t —> — oo and h — > do not commute. Denote U h (t) the 
unitary group associated to the linear semi-classical Schrodinger equation (/ = 
in the above equation). Then for fixed h > 0, 

Fix h in the above asymptotics. If / has moderate growth as in Remark ll.2l fcubic 
nonlinearity in space dimension two or three for instance), then there is scattering 
for 16.2J1 with h fixed, and 

for some concentrating profile where the powers of h stem from the asymptotics 
for the free operator U h . We saw in the introduction that at least when the non- 
linearity is homogeneous, the conservation of energy rules out such a possibility, in 
the limit h — » 0. 
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It is probably more interesting to try to match with the results of Proposition ^. II 



Comparing (|4.1|l and l|6.1(l . Proposition 14. 21 yields . for -J- < t < =tj 



t n — — T-\-t R ' 

2 

Note however that the term in the phase does not belong to any Sobolev 
space; we would have to adapt the statement of Proposition 14.21 before claiming 
that we have ft independent estimates. We can then try to use Grenier's ideas 
again to extend the lifespan of a h and <& ri to [ ~\ h , T], with suitable ft independent 
estimates. We shall not pursue this point of view here. 

We conclude this section by listing a series of questions that remain: 
Do we have a WKB like description of u £ for some time £ > 1? If yes, with one 
or several phases? 

Is there a caustic for ip h l This is not even clear. Indeed, the initial problem (|1.2(l 
contains a data which causes focusing at one point in the linear case. However, we 
saw above that when reaching the boundary layer of size e 1 , phase and amplitude 
of the solution become coupled in such a way that the geometry of the propagation 
is modified. If by any chance there is no caustic for ip , then we might take T 
arbitrarily large, and hope to get a description of u £ for any time. 

Note also that when WKB asymptotics is valid for (p h , then the nonlinear term 
in the conservation of the energy reaches its maximal order of magnitude. This is 
how we found the parameter 7 in the introduction for an homogeneous nonlinearity. 
In the case of Assumptions ITTl things are similar. Denote 

rv 

F{y) = / / ( 7 ? 2 ) vd-V (f] is a real variable). 
Jo 

Then the generalization of the conservation of energy in l|1.3fl is: 

\\\eV x u £ (t)\\ 2 L 2+E- k I F (e k/2 \u £ (t,x)\) dx = const. = 0(1) ~ ^||xa |||2 . 



The nonlinear term in the above energy is exactly J F (\ip h (^r, 2;) |) dx, and is of 
order 0(1) for, say, £ [^?, ^f], where WKB asymptotics for ip h stems from 
Proposition 14. II (For — > —00, dispersive properties of ip h make the nonlinear 
term small.) On the other hand, we saw that if there is a caustic for (f h , then it 
cannot be reduced to a (single) point. 

Finally, the apparent instability discussed in Section suggests that computing 
reliable numerical simulations to understand the asymptotic behavior of u £ is a 
challenging problem, even before a caustic is formed, if there is any. Understanding 
the behavior of (a h ,4> h ) and (a fi, ,<E> fi ), rather than working on u £ directly, would 
certainly be more reasonable. 
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